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' Using the quark- meson coupling (QMC) model, we study the longitudinal re- 

O ! 

' spouse function for quasielastic electron scattering from nuclear matter. In QMC 



the coupling constant between the scalar (o") meson and the nucleon is expected to 
decrease with increasing nuclear density, because of the self-consistent modification 
of the structure of the nucleon. The reduction of the coupling constant then leads to 
a smaller contribution from relativistic RPA than in the Walecka model. However, 
since the electromagnetic form factors of the in-medium nucleon are modified at the 
same time, the longitudinal response function and the Coulomb sum are reduced by 
a total of about 20% in comparison with the Hartree contribution. We find that the 
relativistic RPA and the nucleon structure variation both contribute about fifty-fifty 
to the reduction of the longitudinal response. 
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There is still considerable interest in the longitudinal response for quasielastic electron 
scattering. Within the framework of nonrelativistic nuclear models and the impulse ap- 
proximation, it is very difficult to reproduce the observed, quenched longitudinal response 
functions In the mid '80s, several groups calculated the longitudinal response function 
using Quantum Hadrodynamics (QHD) (i.e., the Walecka model). They argued that 
the contribution of the relativistic random phase approximation (RRPA), which includes 
vacuum polarization, is very important in reducing the Coulomb sum rule i, i below the 
sum of the squares of the nucleon charges in the nucleus. There have also been several 
other attempts to study the longitudinal response in nonrelativistic approaches 0. 

On the other hand, the nucleon has internal structure, and it is nowadays expected 
that this structure should be modified in a nuclear environment ^ . This is closely related 
to the issue of chiral restoration in QCD. In QHD nuclear matter consists of point-like 
nucleons interacting through the exchange of point-like scalar (a) and vector (u) mesons. 
It would clearly be very interesting to investigate the quenching of the longitudinal re- 
sponse function in a relativistic framework, including in addition, the structural changes 
of the nucleon in-medium. 

Recently, we have developed a relativistic quark model for nuclear matter, namely, the 
quark- meson coupling (QMC) model 0], which could be viewed as an extension of QHD. 
However, in QMC the mesons couple to confined quarks (not to point-like nucleons) and 
the nucleon is described by the MIT bag model. This model yields an effective Lagrangian 
for a nuclear system , which has the same form as that in QHD with a density dependent 
coupling constant between the a and the nucleon (N) - instead of a fixed value. Indeed, 
from the point of view of the energy of a nuclear system, the key difference between QHD 
and QMC lies in the a-N coupling constant, Qs- Although this difference may seem subtle, 
it leads to many attractive results [0, We have already applied this model to various 
nuclear problems |p. Here we use it to study the effect of nucleon structure variation in 
the longitudinal response function from nuclear matter. 

First, let us briefly review the calculation of the longitudinal response function for 
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quasielastic electron scattering from (iso-symmetric) nuclear matter in QHD. The starting 
point is the lowest order polarization insertion, H^u, for the u meson. This describes 
the coupling of a virtual vector meson or photon, of momentum q, to a particle-hole or 
nucleon-antinucleon excitation: 

= -Wl J ^jT[G{kh,G{k + g)7.], , (1) 

where G{k) is the self-consistent nucleon propagator (with momentum k) in relativistic 
Hartree approximation (RHA) given as 
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Here fc*^ = - g^V'^ , k) is the mean value of the u field). El = \J k'^ + M*^ (M* is 
the effective nucleon mass in matter) and fc^? is the Fermi momentum. Using the nucleon 
propagator we can separate the polarization insertion into two pieces: one is the density 
dependent part, 11^^, which involves at least one power of G/j, and the other is the vacuum 
polarization insertion, 11^^,, which involves only Gp- The former is finite, but the latter is 
divergent and must be renormalized. We choose to renormalize such that Ii^j^{q) vanishes 
at = and M* = M (where and M are respectively the free masses of the u 



meson and the nucleon). We then find |]10 
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We assume that the isospin degeneracy of the vacuum is 2. For n£,, the exphcit, analytical 



expressions can be found in Ref. (also see Ref. |[T0[|). 

In the Hartree approximation, where only the lowest one nucleon ring is considered, 
the longitudinal response function, , measured in electron scattering is simply given 
by 

„ ) = - {^!Mf\ 9„n,(,). (6) 

Here Z is the nuclear charge, pb the nuclear density, IIl{= ^33 — Hoo) the longitudinal 
component of the polarization insertion (we choose the direction of q as the 2;-axis) and 
GpE is the proton electric form factor, which is usually parametrized by a dipole form in 
free space: 



1 

[l + Q~/0.7iy 



with the space-like momentum transfer, = —qf^, in units of GeV^. For this initial 
investigation we omit a small (and rather complicated) contribution from the anomalous 
moments 0, in order to concentrate on the role of the variation of the structure of 
the nucleon. Since the vacuum polarization is real in the space-like region there is no 
modification of the Hartree response from this term. 

The RRPA for the longitudinal component of the polarization insertion, Ilf;^^, involves 
the sum of the ring diagrams to all orders. This summation has been discussed by many 
authors [0, ^, 0, [Tl|, . It involves a-uj mixing in the nuclear medium, and is given by 



nr^(g) = [(1 - AoH,)H^ + AoUU/e,, (8) 
where e^, is the longitudinal dielectric function 

ez. = (1 - doUL){l - AoH,) - {ql/q')AodoUl, (9) 

with q = \q\, and the free meson propagators for the a and u mesons are respectively 

1 1 

Ao(g) = — 5-— and do{q) = — 5——, (10) 
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where rricr is the a meson mass. Here lis and 11^ are respectively the scalar and the time 
component of the mixed polarization insertions: 
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Ti[G{k)G{k + q)], 
Tr[G(A;)7°G(fc + g)]. 



(11) 

(12) 



The scalar polarization insertion can be again separated into two pieces. The density 
dependent part is finite and the explicit expression can be found in Ref. |n|. Because it 
does not involve Go, the vacuum component, Ilf , is, of course, divergent and once again 
we need to renormalize it. First, we introduce the usual counter terms to the Lagrangian, 
which includes terms quadratic, cubic and quartic in the a field, as well as wavefunction 
renormalization p[. To get the "physical" properties of the a meson in free space, we 
impose the following condition [|10[] : 



nf(g;,M* = M) 

Then, we find 
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where Zs = 1 — For the mixed polarization insertion there is no vacuum polarization 
and it vanishes at zero density. (The explicit form can be also found in Ref. |1TT| .) 

As QHD involves only isoscalar mesons, the isovector RRPA response is the same 
as the Hartree one, eq.(|^). This implies that the vacuum polarization only affects the 
isoscalar response. It remains to study the effect of isovector mesons. (In the isovector 
part the rho meson coupling (without vacuum polarization) was studied in Ref. It 
reduces Sl slightly.) Since the longitudinal response is half isoscalar and half isovector. 



the longitudinal response function in RRPA is given by 

(15) 
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Several authors |] have calculated the longitudinal response function using this RRPA 
polarization, and reported that it is very important in reproducing the observed experi- 
mental data [0. 

Now we are in a position to discuss the effect of changes in the internal structure of 
the nucleon in-medium. In order to do so, we consider the following modifications to the 
QHD approach: 

(1) meson- nucleon vertex form factor 

In QHD the interactions between the mesons and nucleon are point-like. However, 
since both the mesons and nucleon are composite they have finite sizes. In the region of 
space-like momentum transfer the finite-size effect reduces the meson-N coupling. As the 



simplest example, one could take a monopole form factor ||14| at each vertex: 



with a cut off parameter Ajy = 1-5 GeV. In principle, one could self-consistently calculate 
the form factor within QMC However, as such changes are not expected to make a big 
difference, we use eq.(|16|) in the following calculation. 
(2) modification of the proton electric form factor 

Recently we have studied the electromagnetic form factors of the nucleon, not only 



in free space |]T5[ but also in a nuclear medium, using the QMC model [0 (see also 



Ref. IT^l)- Because the confined quark feels an attractive force due to the a, the quark 
wave function is modified in a nuclear medium. The ratio of the electric form factor of 
the proton in medium to that in free space, Gpe{pb,Q^)/Gpe{Q'^), is shown in Fig. 3 in 
Ref. |]TB|. From the figure we can see that the ratio decreases very linearly as a function 
of Q^, and that it is accurately parametrized at ps = Po (= 0.15 fm~'^, the normal nuclear 



matter density) as 



RpEiPo, Q') = - 1 - 0.26 X Q\ (17) 



This implies that the (electric) rms radius of the proton at po swells by about 5.5 % (for 
more details, see Ref. |1^). Since the bag model reproduces the form factor measured in 
free space very well and the latter is well described by eq.(|^), the in- medium proton 
form factor can be represented as G'ps(Q^) x Rpe{pb,Q'^^- 
(3) density dependence of the coupling constants 

In QMC the confined quark in the nucleon couples to the a field which gives rise to 
an attractive force. As a result the quark becomes more relativistic in a nuclear medium 
than in free space. This implies that the small component of the quark wave function, 
ipq, is enhanced in medium ||^, ^. The coupling between the a and nucleon is therefore 
expected to be reduced at finite density because it is given in terms of the quark scalar 
charge, J^^g dVipgipg . On the other hand, the coupling between the vector meson and 
nucleon remains constant, because it is related to the baryon number, which is conserved. 

To study the longitudinal response of nuclear matter, we first have to solve the nuclear 
ground state within RHA. In QHD the total energy density for nuclear matter is written 



as 



IC 
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where So has the usual form (in RHA), given in Ref. 0. Note that in Ref. the 
renormalization condition on the nucleon loops is imposed at qfi=0. The second term on 
the r.h.s. of eq.(|l8D 0] occurs because we chose the renormalization condition for the a 
at = ml (see eq.(|l3D). As measureable quantities cannot depend on this choice, our 
model gives the same physical quantities as those of Ref. . 

To take into account the modifications (1) and (3), we replace the cr- and uj-N coupling 
constants in eq.(|T^) by 

gs ^ ^7s(Pb) X F^(Q2), (19) 
^ g.xF^iQ^), (20) 



^ The nucleon-antinuclcon excitation in H^^ contributes to the photon-nuclcon vertex as a RRPA 
correction, which may in principle lead to double counting for the form factor. However, we ignore this 
correction because it is very small fl . 



7 



where the density dependence of gsips) is given by solving the nuclear matter problem 
self-consistently, using the MIT bag for the nucleon model (see Ref. As in QHD, 

we have two adjustable parameters in the present calculation: gs{0) (the cx-N coupling 
constant at = 0) and g^j. 

Requiring the usual saturation condition for nuclear matter, namely S/pb — M = 
—15.7 MeV at po, we determine the coupling constants gs{0) and g^ (fi'KO) = 61.85 and 
g^ = 62.61). In the calculation we fix the quark mass to be 5 MeV, rricr = 550 MeV 
and rrii^ = 783 MeV, while the bag parameters are chosen so as to reproduce the free 
nucleon mass (M = 939 MeV) with the bag radius Rq = 0.8 fm (i.e., B^/^ = 170.0 
MeV and z = 3.295 ^ g). This yields the effective nucleon mass M*/M = 0.81 at po 
and the incompressibility K = 281 MeV. (We do not consider the possibility of medium 
modification of the meson properties [Q in the present work.) 

1 

0.95 
0.9 

o 

'■B 0.85 
DC 

0.8 
0.75 
0.7 

0.5 1 1.5 2 

Figure 1: Density dependences of gs{pB) / gs{^) and M* /M. The solid curve is for the ratio 
of the coupling constants, while the dotted curve is for the ratio of the nucleon masses. 

Now we present our main results. First, in Fig. |l|, we show the density dependence of 
the coupling constant. At po, gs decreases by about 9%. The effective nucleon mass is 
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Figure 2: Longitudinal response functions in QMC with rriq = 5 MeV. We fix g = 550 
MeV and pb = Po- The dotted curve is the result of the Hartree approximation (see 
eq.(P)), where the effective nucleon mass is given by QMC and the proton electric form 
factor is the same as in free space. The dashed curve is the result of the full RRPA, 
without the modifications (1) and (2) (i.e. = 1 and RpE = 1). The dot-dashed curve 
shows the result of the full RRPA with the meson-N form factor but RpE = 1. The upper 
(lower) solid curve shows the result of the full RRPA for ruq = 5 (300) MeV, including all 
modifications. 

also shown in the figure. 

Next, we show the longitudinal response function in QMC. Using the density de- 
pendent coupling constant, the meson-N form factors (see eqs.(|19|) and (pOD) and the 
in-medium proton electric form factor, we can calculate the longitudinal response of nu- 
clear matter. The result is shown in Fig. |^. Because of the density dependent coupling, 
9s{pb), the reduction of the response function due to the full RRPA (the dashed curve in 
the figure) from the Hartree result (the dotted curve) becomes much smaller than that 
in QHD. On the other hand, the modification of the proton electric form factor is very 
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Figure 3: Coulomb sum, C{q)/Z, at po in the Hartree approximation with RpE = 1 (the 
dotted curve) or the full RRPA with all modifications (the solid and dot-dashed curves 
are for = 5 and 300 MeV, respectively). 

significant, yielding a much bigger reduction in the response (see the upper solid curve). 
We can see that the effect of the meson-N form factor enhances the longitudinal response 
(see the dot-dashed curve), but it is not large. 

It is also interesting to see the quark mass dependence of the longitudinal response. 
As an example, we consider the case of rrig = 300 MeV, which is a typical constituent 
quark mass. For = 300 MeV and Rq = 0.8 fm, the coupling constants required to 
fit the saturation properties of nuclear matter are: gg{0) = 68.69 and g"^ = 84.24, and 
the effective nucleon mass at po and the incompressibility become 723 and 345 MeV, 
respectively. Using these parameters we show the result for the longitudinal response 
(the lower solid curve) in Fig. |^. In comparison with the case rriq = 5 MeV, it is a little 
smaller and the peak position is shifted to the higher energy transfer side. This may be 
due to the smaller effective nucleon mass in the case mg = 300 MeV than when = 5 
MeV. 
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The integrated strength of the longitudinal response (or the Coulomb sum), C{q), 



C{q)= / dqoSL{q,qo), (21) 
Jo 

is shown in Fig. ^ as a function of three-momentum transfer, q. For high q, the strength 
is about 20% lower in the full calculation than for the Hartree response. For low q, the 
full calculation with the constituent quark mass remains much lower than the Hartree 
result, while in case of the light quark mass it gradually approaches the Hartree one. This 
difference is caused by that the effective nucleon mass for rUg = 5 MeV being larger in 
matter than that for rrig = 300 MeV. The 20% reduction found here is a little smaller 
than the value of approximately 30% found in QHD |^, |[. 

We would like to emphasize that these calculations are for nuclear matter and can- 
not be directly compared with the experimental data. Furthermore, there still remain 



discrepancies and uncertainties in the present experimental results |18, 19 



Finally, we comment on the transverse response from nuclear matter. In Ref. |16] we 
can see that in QMC the modification of the nucleon magnetic form factor in-medium 
is very small: the calculated decrease in the proton (neutron) magnetic form factor is 
about 1.5% (0.9%) at po- Therefore, one would expect the total change in the transverse 
response caused by RRPA correlations and the effect of the variation of the structure of 
the nucleon to be much smaller than in the longitudinal response. This is certainly what 
one needs in order to fit the experimental data 

In summary, we have calculated the longitudinal response of nuclear matter using the 
QMC model. The reduction of the a-N coupling constant with density decreases the 
contribution of the RRPA, while the modification of the proton electric form factor in 
medium reduces the longitudinal response considerably. The longitudinal response, or 
the Coulomb sum, is reduced by about 20% in total, with RRPA correlations and the 
variation of the in-medium nucleon structure contributing about fifty-fifty. In the near 
future we hope to extend this work to calculate the longitudinal and transverse response 
functions for finite nuclei, using local density approximation, in order to compare our 
results directly with new experimental results |]T9 . 
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